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^ . ABSTRACT 

f ^ ' Using photometric redshifts we determine the galaxy population of the clusters of 

^ ! galaxies C10016+16 at z = 0.55, C11600+41 at z = 0.54, C11601+42 at z = 0.54 

■ and MS1008~1224 at z = 0.31. Comparing the clusters, we find no evidence for a 

\ universal shape of the total luminosity function (LF) at these redshifts. When dividing 

. the LFs into spectral types, we find that the LF of the early- type galaxies alone can 

be described by a Gaussian, while the LF of the late-type galaxies is well fitted by 
a Schechter function, suggesting that the separate LFs for different populations may 
be universal. The difference in the total LFs can mainly be attributed to the varying 
relative normalisation of these populations, implying that clusters with an abundant 
population of late- type galaxies also have steeper faint-end slopes. In MS1008-1224 we 
detect a faint blue population that dominates over a population with colours consistent 
I with dwarf ellipticals, opposite to clusters at lower redshift. 

, Compared to low-redshift clusters, we find that a general fading of the late-type 

C\| ' population by 2 mag and the early-type population by 1 mag describes the 

\ evolution from z = 0.55 to z = well. 

. As a consequence of the different early-type and late-type LFs and their depen- 

' dence on cluster radius, the fraction of IdIuc cluster galaxies, as measured by the 

fH , Butcher-Oemler effect, differs between the clusters and depends on limiting magni- 

i' tude and radius. 

We find a correlation between the dwarf-to-giant ratio and the surface density, 
indicating that the high density environment in the cluster cores is hostile to dwarf 
galaxies. 

' Key words: galaxies: clusters: general - galaxies: distances and redshifts - galaxies: 

luminosity function, mass function - cosmology: observations. 

1 INTRODUCTION higher accretion rate of star forming field galaxies at higher 

,rT-,io^i/ r redshift (Bower 1991), possibly combined with a general in- 

Ever smce the results of Butcher & Oemler (1984, hereafter . . ^ , . ■ ri , , ■ ^ i • i /t^- r • 

^ crease m star formation m field galaxies at niglier z (Diaferio 

B084), showing that the fraction of blue galaxies in clusters , , . rj_i ^-.t-mm j i - j_ i ■ i 

-ji T n J n c • u u et al. 2001). A consequence of the CDM model is that high-z 



o 



increases rapidly between z = Q and z ~ 0.5, it has been 
clear that there is strong evolution in cluster galaxy popu- 



clusters assemble during a much shorter time interval than 

comparably rich clusters today (Kauffmann 1995). This gen- 
lations with redshift. These results have been confirmed by , . . , ^ i i j_i i i ■ r 

eral scenario is also supported by the observed increase of 
subsequent surveys (e.g. Rakos & Schombert 1995: Ellmg- • i • i i i i ^ / i i ^ i 

t'^i ot, ^, , mergers m nign-redslnft clusters (van Dokkum et al. 1999). 

son et al. 2001; Kodama & Bower 2001; Margonmer et al. j_i at i i ■ i j_- i t~, i • i • ■,nnc\ i 

° lurtnermore, JN-body simulations (e.g., Dubmski 1998) snow 

2001). A general picture explaining the over-all features of ^ i ■ , ^ ^ i ^ i • i-nr^r^\ j_ n r 

, ' , , that the brightest cluster galaxies (BOG) naturally torm via 

cluster formation and evolution is given withm the cold dark . , . , • , 

,r^rs^,\ ,. ,. , n • . / t, mergBrs lu a hierarchical scBuario. 

matter (CDM) hierarchical clustering scenario (e.g. Bower 

1991; Kauffmann & White 1993; Kauffmann 1995; Baugh, 

Cole & Frenk 1996). A large number of observations and mi i-rr • i- i i- ^ -i j_- i ^ i j_ 

' ° ihe difference m radial distribution between early-type 

N-body simulations now support this scenario. The blue- u i-ju^i 

^'^ and late-type galaxies m clusters can be explained by the 

ing with redshift of the cluster galaxies is explained by a , . i • i j i -r • r n- c u i • i ■ i 
° o 1- J hierarchical model il mtallmg field galaxies, which are pre- 

dominantly blue star-forming galaxies, have their star for- 
mation truncated as they fall deeper into the cluster poten- 
* E-mail: tomasQastro.su.se. Present address: Space Telescope tial (Balogh, Navarro & Morris 2000; Diaferio et al. 2001). 
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ment also aflect the star formation, transforming the spec- 
tral type of the galaxies. 

One of the most important probes of cluster evolution is 
the luminosity function (hereafter LF), describing the num- 
ber of galaxies per magnitude bin in a cluster. Binggeli, 
Sandage & Tammann (1988) showed that the total LF for 
the Virgo cluster rises as a Gaussian at bright magnitudes 
and then flattens to a "plateau" (or even a slight decrease), 
before it turns steep at faint magnitudes. LFs with similar 
shapes arc found in other nearby clusters, e.g. Coma (Tren- 
tham 1998a). The existence of a numerous population of 
faint dwarf galaxies in clusters, more common than in the 
field, is verified by a number of observations (e.g.. Smith, 
Driver & Phillipps 1997; Wilson et al. 1997; Phillipps et al. 
1998; Trentham 1998a; Yagi et al. 2002). There are, how- 
ever, also clusters that do not show this steep increase at 
faint magnitudes (Trentham 1998b, 1998c). Most likely, this 
is a consequence of different evolution depending on e.g., 
richness and epoch of formation. 

Conselice, Gallagher & Wyso (2001) discuss different 
scenarios that could lead to the dwarf population in local 
clusters. By comparing models with the observed kinematic 
and spatial properties of a number of Virgo cluster galaxies, 
they suggest that a transformation of spirals into dwarf ellip- 
ticals (dEs) by "galaxy harassment" (Moore, Lalce & Katz 
1998) is occurring. Some dEs may, however, have formed 
outside the cluster in galaxy groups, which were later ac- 
creted. These dEs may therefore be as old as the cluster 
ellipticals. The fading of these cluster galaxies into todays 
dwarf population is discussed by Wilson et al. (1997) and 
Small et al. (1998). 

Besides measuring the total LF, it is of high interest 
to investigate the LFs of different galaxy populations and 
their evolution with redshift. While previous studies of type- 
dependent LFs have focused on morphological types, e.g., 
Binggeli et al. (1988) and Andreon (1998), here we divide 
galaxies into spectral types characterised by their colours, 
as described in Section 3.3. In this paper, when we refer to 
galaxies as early-types and late-types we mean their spectral 
type, while morphological types are denoted by their type 
in the Hubble sequence, e.g., ellipticals and spirals. 

A consequence of the hierarchical scenario is that clus- 
ters observed at high redshift {z ^ 0.5), where the infall of 
field galaxies is assumed to pealc (Bower, Kodama & Ter- 
levich 1998), should have a higher fraction of late- type galax- 
ies compared to local clusters. This should lead to a differen- 
tial evolution between the late-type and early-type cluster 
LFs with redshift, with the late- type LF shifted towards 
brighter magnitudes at increasing redshift, as compared to 
the early-types. 

Because it is observationally costly, if at all possi- 
ble, to obtain redshifts for the numerous faint population, 
most observational efforts to study galaxy populations in 
high-0 clusters have concentrated on the brighter galaxies 
{Mb ^ —19). Recently, however, photometric redshifts have 
been demonstrated to be a powerful tool for studying espe- 
cially the faint population of cluster galaxies. Even though 
the accuracy of photometric redshifts can not be compared 
to spectroscopic redshifts, they can provide a reliable deter- 
mination of clusters membership, and therefore significantly 
reduce the necessary amount of field galaxy subtraction. In 
a previous paper (Dahlen, Fransson & Naslund 2002, here- 



after DFN02), wc demonstrated the use of photometric red- 
shifts in a study of the population of the intermediate rich 
cluster C11601-I-42 aX z = 0.54. The photometric redshift 
selection minimizes the amount of background subtraction 
needed. Internal properties of the cluster, such as radial dis- 
tribution and luminosity function, both for the total cluster 
population, as well as for different populations separately, 
can therefore be determined. Here we extend our previous 
study of C11601-f42 to include two additional rich clusters, 
C10016+16 at z = 0.55 and MS1008-1224 at 2 = 0.31. We 
also include observations of a fourth poor cluster, C11600-I-41 
at z = 0.54. 

Throughout this paper wc assume a Hubble constant 
Ho = 50 km s^^ Mpc^^, and a cosmology with Om = 0.3 
and SIa = 0.7, unless anything else is assumed. Magnitudes 
are given in the Vega based system. 



2 THE DATA 

2.1 Observations 

The positions of the four clusters and a blank field, used for 
background subtraction, are listed in Table 1. All observa- 
tions, except those of MS1008 1224, wore carried out with 
the 2.56-m Nordic Optical Telescope (NOT) and the An- 
dalucia Faint Object Spectrograph and Camera (ALFOSC) 
during six observing runs between 1997 and 2001. The clus- 
ters and the background field were observed in four filters, B, 
V, R and /, and additional observations in the U filter were 
obtained for C11601-I-42. Observations were performed un- 
der photometric conditions. The seeing in the images varies 
between 0".70 and 1".15. A complete log of the observations 
is given in Table 2. 

The observations of MS1008-1224 were carried out by 
the Science Verification Team at ESO using FORS at the 
VLT. The cluster was observed in B, V, R, I FORS Bessel 
filters with exposure times ~ 1 — 1.5h. The images were re- 
duced using the iraf package. The seeing in the final coad- 
ded frames is (0".72, 0".65, 0".64, 0".55) in the {B,V,R,I) 
filter. Further details can be found at the ESO web site^. 
Information on redshifts, galactic extinction, observed area. 
X-ray luminosities and velocity dispersion of the clusters is 
given in Table 3. 

2.2 Data reductions 

The data obtained with the NOT were reduced using the 
IRAF package. Bias subtraction and flat-fielding were made 
in a standard manner. For the /-band we constructed a 
fringe- frame after removing objects from the science images. 
The fringe-frame was then subtracted from each science im- 
age, scaled to the appropriate background level. The images 
were corrected for atmospheric extinction, aligned, and fi- 
nally combined. Calibration was done using standard stars 
from Landolt (1992). The galaxies were corrected for galac- 
tic extinction according to Schlegel, Finkbeiner & Davis 
(1998). A description of data reduction and calibration of 
the MS1008-1224 images is given at the ESO web page. 

t URL http://www.hq.eso.org/science/utlsv 
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Table 1. Positions of the observed clusters and back ground field given in J2000 coordinates. 

Object RA Dec 

C10016+16 00''18'"33^.3 16° 26' 36" 

C11600+41 16''02'"06M 41° 01' 23" 

C11601+42 16''03'"09*.8 42° 45' 18" 

MS1008-1224 10''10'"34«.1 -12° 39' 48" 

Blank field 16''08'" 54^.0 41° 34' 00" 



Table 2. Log of observations. 



Object 


Filter 


Obs. date 


Exp. time 


# of exp. 


Seeing 


liT (mag arc 


C10016+16 


B 


Aug 00 


18000s 


20 


0".92 


27.7 




V 


Aug 00 


17100s 


19 


1".03 


27.3 




R 


Aug 00 


7200s 


8 


0".73 


26.6 




I 


Aug 00 


16200s 


36 


0".74 


26.0 


C11600+41 


B 


Jun 01 


18000s 


20 


0".89 


27.8 




V 


Jun 01 


12600s 


14 


0".86 


27.2 




R 


Jun 01 


8100s 


9 


0".89 


26.7 




I 


Jun 01 


18000s 


30 


0".83 


26.2 


C11601+42 


U 


Apr/ Jun 99 


25200s 


14 


0".89 


26.7 




B 


Jun 98 


12600s 


14 


0".84 


27.5 




V 


Apr 99 


16200s 


18 


1".02 


27.3 




R 


Jun 97 


7200s 


8 


0".75 


26.8 




I 


Jun 98/ Jun 99 


14400s 


24 


0".76 


26.2 


Blank field 


B 


Aug 00 


16200s 


18 


1".15 


27.6 




V 


Jun 99 


16200s 


18 


1".02 


27.0 




R 


Jun 97 


7200s 


8 


0".70 


26.8 




I 


Aug 00 


16200s 


36 


0".73 


26.8 



Photometry was obtained using the FOCAS package 
(Jarvis & Tyson 1981; Valdes 1982; Valdes 1993). The de- 
tection limit was set to Sct of the sky noise, and a minimum 
detection area corresponding to the seeing-disc was used. For 
each object and filter, we made a catalogue listing isophotal 
magnitude, aperture magnitude, position and area. When 
calculating the aperture magnitude we first smoothed the 
observations to the seeing of the filter with the worst seeing 
of each cluster, and then used a constant aperture size with 
a radius corresponding to this seeing. Finally, a combined 
catalogue for each object was made by using the positions 
in the i?-catalogue and matching these with the positions 
in the other catalogues. 

The completeness was tested with simulations where 
we added artificial galaxies with different magnitudes and 
radial profiles to the real images. Using the same detection 
procedures as for the real data we find that 100 per cent 
of the galaxies are detected down to ruR = 25 for all 
our objects. Using simulations we also calculate corrections 
that should be applied when calculating total magnitudes 
from the observed isophotal magnitudes. These simulations 
are described in Naslund, Pransson & Huldtgren (2000) and 
DFN02. 

In Table 4 we list the number of objects in each cata- 
logue to a limiting magnitude mn < 25. We also list the 
number of objects having aperture photometry in at least 
four filters (i.c four or five filters for C11601+42, and four 
filters for the other objects), three filters, as well as in two 
or one filter only. Objects identified as stars by visual in- 
spection of the psf, are excluded. 



3 PHOTOMETRIC REDSHIFTS 

An extended discussion on the use of photometric redshifts 
applied to clusters of galaxies is given in DFN02. Here we 
present a brief summary of the technique. For every object 
we minimize the expression 

2,. ^ ^ {nii - {Ti{t, z) + mc)f 
X {t,z,mc) = 2_^ ^ (1) 

i * 

where rm and ffi are the observed magnitudes and uncer- 
tainties in filter i, respectively. Ti{t,z) is the magnitude in 
the i-filter of template t, redshifted to z. This magnitude is 
calculated by convolving the filter curve and the quantum ef- 
ficiency of the detector with the galaxy template. The quan- 
tity ma is a constant, which fits the apparent magnitude of 
the template galaxy. 

A set of ten different templates are used. We construct 
these by interpolations between the four observed galaxy 
templates given by Coleman, Wu & Woodman (1980). These 
represent E, Sbc, Scd and Im galaxies. Absorption due to 
intergalactic Hi clouds is treated as in Madau (1995). We 
also include eight stellar templates of M-dwarfs taken from 
Gunn & Stryker (1983). 

A modification compared to DFN02 is that we here use 
a "Baycsian" approach (e.g. Kodama, Bell & Bower 1999; 
Benitez 2000) . This method allows the incorporation of pre- 
existing knowledge about the galaxies into the photometric 
redshift determination. This is illustrated in Fig. 1, where 
we in the top panel show the probability distribution do- 
rived from the chi-square fit according to Eq.(l), for a galgixy 
with known spectroscopic redshift z = 0.31. There are ten 
local minima in the chi-square fit, resulting in ten probabil- 
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Table 3. Properties of the cluster sample. 



Cluster 


Redshift 




Observed area 


Diameter^ 


Lr (0.3-3.5 keVO 










s(]. arcuiiii 


Mpr 


10-l-i <.r.t^ s 


km s-1 


C10016+16 


0.546'= 


0.19 


32.5 


3.1 


34.7'' 


ires'* 


C11600+41 


0.540' 


0.04 


28.9 


2.9 


<2.0''-'' 




C11601+42 


0.539'= 


0.03 


30.4 


3.0 


2.1=* 


1166'' 


MS1008-1224 


0.306-'' 


0.23 


32.4 


2.2 


5.89 


1054'' 



Notes: 

a) Schlegel et al. 1998 

6) Cluster diameter covered by observations. 

c) Dressier & Gunn 1992 

d) Small et al. 1997. 

e) Oke, Gunn & Hoessel 1996. 

/) Lewis ct al. 1999. 

g) Gioia & Luppino 1994. 

h) Carlbcrg ct al. 1996. 

i) Henry et al. 1982. 

j) Converted from energy band 0.5-4.5 keV assuming electron temperature 4 keV 



Table 4. Number of objects with nijj < 25 in the different images, and the number of objects with aperture 
photometry in four to five filters, three filters, and one or two filters only. Stars are excluded. 

Image No. of Objects > 4 filters 3 filters 1-2 filters 



C10016+16 1430 1154 172 104 

C11600+41 809 791 13 5 

C11601-I-42 1199 1034 106 59 

MS1008-1224 1665 1544 104 17 

Fieldl 945 781 87 77 



ity maxima, divided into two groups, one at low, and one at 
high redshifts. The best-fitting template is in this case an 
Sa galaxy with redshift z = 0.20, but almost as good fits 
are achieved for the three other templates with z < 0.5, 
representing E to Sbc galaxies. There are also maxima with 
high probability at a ~ 2.8 — 3.5. The degeneracy of the 
high- and a low-redshift groups arise because the Lyman- 
break falls approximately between the same filters at « ~ 
3.5, as the 4000-A break does at ^: ~ 0.3. 

The absolute magnitude for this galaxy at z ~ 3.5 
would, however, be unrealistic, Mb ~ — 28, whereas 
at z ~ 0.1 and z ~ 0.3 we get Mb ~ — 18 and 
Mb — 20.5, respectively. To account for this one could 
weight the probability function with an expected LF, an ap- 
proach used by Kodama et al. (1999). However, since we are 
here interested in determining the LF, we can not use this 
method. Instead, we use an exponential cut-off two magni- 
tudes brighter than Mb, whore Mb is determined by a fit 
to the Schechter function without any weighting. The effect 
of introducing this cut-off at bright magnitudes is mainly to 
suppress the false peaks introduced by the misidentification 
between the Lyman-break and the 4000-A break. An alter- 
native would bo to truncate the procedure at e.g. z = 2. 
However, at ttir ~ 25, galaxies at « ^ 2 with "normal" ab- 
solute magnitudes are expected; a cut-off could here lead to 
a misidentification of these galaxies. Furthermore, the vol- 
ume element increases with redshift (up to z ^ 1.8, dotted 
line in the top panel of Fig. 1). Therefore, of two maxima 
with the same probability from the chi-square fit, the higher 
redshift should for constant comoving density, be more likely 
due to the larger volume. 



The lower panel of Fig. 1 shows the probability distri- 
bution for the same galaxy after multiplying the result from 
the top panel with the probability distribution derived after 
applying an exponential cut-off at bright magnitudes and 
a weighting by the volume element. The distributions are 
normalised to unit area. 

The peaks at 2; > 2.5 have now disappeared, and the 
lowest redshift peaks are suppressed due to the smaller vol- 
ume element. The maximum probability is now for a Sb 
galaxy at z = 0.31, which matches the spectroscopic red- 
shift. 

The "Bayesian" approach we use mostly affects the de- 
termination of redshifts for the cluster MS1008-1224 at z 
= 0.31. Liu & Green (1998) show that when using photo- 
metric redshifts there is a risk of misidentification between 
Sbc galaxies at 2: ~ 0.05 and Scd/Im galaxies at z ^ 0.3. 
From our example above we also note that ~ Sb galaxies at 
z ~ 0.3 can be misidentified as E galaxies at z ^ 0.1. 

When estimating the dispersion between spectroscopic 
and photometric redshifts for 61 galaxies in the field of 
MS1008-1224, we find cr^ ~ 0.051, when we use the 
"Bayesian" approach, compared to az ~ 0.11 without 
it. The increased dispersion is mainly duo to a few galax- 
ies achieving large errors. Excluding six galaxies reduces the 
dispersion to az ~ 0.058. For the clusters at z ^ 0.55 the 
effect of including this weighting is marginal. Here we get 
Uz ^ 0.054 with the "Bayesian" approach and tj^ ^ 0.060 
without, when calculating the dispersion between the pho- 
tometric and spectroscopic redshifts for 38 galaxies in the 
field containing C10016-I-16. 
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Figure 1. Chi-squcire probability function for the template fitting 
method, for a galaxy with known spectroscopic redshift z = 0.31. 
At each redshift the probability of the galaxy template that has 
the minimum chi-square, i.e the highest probability is plotted. 
The top panel shows the probability function without considering 
absolute magnitudes or volume clement. The dotted line shows 
the redshift dependence of the volume element. In the bottom 
panel we use a cut-off at bright magnitudes, and take the volume 
element into account. For some of the peaks we plot the associated 
galaxy type. The distributions are normalised to unit area. 



The example shown in Fig. 1 was picked to illustrate 
the "Bayesian" method, and has an unusual high number of 
probability peaks with similar strengths. Most galaxies have 
a more clearly defined primary peak. 

Possible systematic effects when using photometric red- 
shifts are discussed in DFN02. In summary, we find a 
marginal increase in the photometric redshift errors for the 
faintest galaxies (m_R^ 24) due to increased photometric er- 
rors. These errors mostly affect late type galaxies. In num- 
bers, we find a possible increase in redshift uncertainty by 
Aaz ~ 0.01 due to systematic effects at mij^ 24. This may 
lead to an underestimate of the number of cluster galaxies in 
the faintest bins by ~5%. This result shows that systematic 
errors are unlikely to dominate over statistical errors. 



3.2 Comment on background counts 

To estimate the background contamination we use BVRI 
photometry of our blank field. In DFN02 we used UBVRI 
photometry of the ESO Imaging Survey (EIS) deep field 
published by da Costa et al. (1998). Number counts to Ji = 
25 for Fieldl and the EIS field yield 983 ± 133 and 845 ± 
120 galaxies, respectively, where the fields are normalised to 
tlic same size as the image of C10016-I-16 (32.5 sq. arcmin). 
Errors represent la and include Poisson statistics and field- 
to-field variance (see DFN02 for a discussion). The resulting 
counts are therefore consistent within the errors. 

When determining the background counts within the 
cluster redshift range, we find a larger deviation between 
the two fields. In the redshift range zciooie ± l-5cr^, we 
find 121 ± 20 galaxies for Fieldl, and 166 ± 26 galaxies 
for the EIS field to i? = 25. These results are marginally 
within errors, but it is likely that there are systematic effects 
responsible for some of the off-set. In particular, the aperture 
magnitudes are determined differently. For our blank field 
we smooth the images to the same seeing before calculating 
the colours, whereas the aperture magnitudes given for the 
EIS field have a varying seeing in the different bands. This 
can affect the colours, especially for small objects that do 
not cover the whole aperture. 

To check the effect of this, we calculate photometric red- 
shifts for the field containing C10016-I-16 without smoothing 
the images to the same seeing. Comparing with the 38 ob- 
jects that have spectroscopic redshifts we find an rms devi- 
ation az — 0.078, which is considerably larger than what 
we found after smoothing, az = 0.054. Therefore, the pho- 
tometric redshifts calculated for the EIS catalogue are likely 
to be less accurate. 

Finally, systematic errors could be introduced by the 
uncertainty in the zero-point magnitudes and the use of dif- 
ferent software when doing the photometry, i.e the use of 
FOCAS for Fieldl, and SExtractor for the EIS field. 

To minimize the risk of introducing systematic errors, 
we use the photometry from our blank field in this analysis, 
since this is derived in the same way as the photometry of 
the cluster images. The use of a different background field 
does, however, not affect the result in DFN02 on C11601-I-42 
more than marginally, i.e. any differences arc within the 
quoted errors. This illustrates the advantage of the photo- 
metric method for examining high-z clusters; the relatively 
small corrections for background contamination makes the 
results less dependent on the background. 



3.1 Photometric redshift catalogue 

We calculate photometric redshifts for all objects with 
ruR < 25 and aperture photometry in at least three filters, 
which corresponds to 92-99 per cent of the total number of 
objects to this limit (Table 4) . The accuracy of the photo- 
metric redshifts is estimated by comparing our result with 
available spectroscopic redshifts. In Table 5 we list the num- 
ber of available redshifts and the rms deviation between the 
photometric and spectroscopic redshifts, az- The galaxies 
used for calculating the rms have photometry in all filters. 
To estimate the increase in errors for objects only having 
photometry in three filters, we also calculate the rms devia- 
tion after excluding the filter with the largest error for each 
galaxy, ai^\ 



3.3 Cluster membership and galaxy classification 

For selecting cluster members we use a redshift range ± 
1.5az around the mean redshift of the cluster, with values of 
az given in Table 5. To account for the increase in dispersion 
for the galaxies only detected in three filters (Table 3), we 
use the corresponding az^^ for these galaxies. Note, however, 
that the information given by a non-detection (i.e., that the 
magnitude of the galaxy is fainter than the magnitude limit 
of the observations) is taken into account in the template 
fitting procedure. Assuming that the errors are Gaussian, 
the 1.5az cut should include 86 per cent of the actual 
number of cluster galaxies to our limit. Taking into account 
the small fraction of the galaxies observed in three filters 



© 2004 RAS, MNRAS 000, 1-16 



6 Dahlen et al. 



Table 5. N is the Number of galaxies with spectroscopic redshifts used to determine the dispersion between photometric and spectroscopic 
redshifts, az is the dispersion using information from all filters, while is the dispersion using only three filters. Also given is the 
resulting number of clusters galaxies within zci ± l-5c7z for the total areas and inside a common radius of 1 Mpc. We also list the density 
of background galaxies within each redshift range. The limiting magnitude is Mg = —17.7, corresponding to mn ~ 25 at z = 0.54 and 
mR ~ 23.4 at z = 0.31. 



Cluster 


^Cl 


N 






Cluster g 


;alaxics 


Background" 












Total area i? < 1 Mpc 


(Mpc-2) 


C10016+16 


0.546 


38'' 


0.054 


0.080= 


463 ± 32 


248 ± 19 


12.7±2.1 


C11600+41 


0.540 


d 


0.065'* 


CIC* 


36 ± 24 


26± 12 


14.6±2.4 


C11601+42 


0.539 


78'' 


0.076 


0.12<= 


332 ± 35 


154 ± 18 


18.8±2.9 


MS1008-1224 


0.306 




0.051 


0.067 


221 ± 20 


173 ± 16 


11.6±2.3 


MS1008-1224S 










319 ± 26 


240 ± 21 


20.6±3.5 



Notes: 

a) Measured within the redshift range defining the different clusters, i.e. zci ± l.bcTz- 
h) Dressier et al. 1999. 

c) Excluding two outliers with large errors. 

d) No spectroscopic redshifts available, except for central galaxy. We assume dispersions 
equal to the mean of the dispersions of the two clusters at similar redshift. 

e) Excluding five outliers. 
/) Ycc et al. 1998. 

g) Number of galaxies to Mg = —16.2. 



only that may result in outliers with large errors, wc estimate 
that we include ~ 80 — 85 per cent of the total number of 
cluster galaxies when applying the l.Scr cut. For a discussion 

on completeness and contamination when selecting cluster 
members using photometric redshift, see Brunricr & Lubin 
(2000). 

In this study we divide the galaxies into late-types and 
early-types, based on the spectral type determined by the 
photometric colours of the galaxies. The division is made 
half way between the E and Sbc templates from Coleman 
et al. (1980), approximately corresponding to a rest-frame 
colour B — V = 0.8. The first category consists of early-type 
galaxies with red colours, which are mostly ellipticals and 
lenticulars, but also include passive spirals with red colours. 
The second category includes late-type spirals and irregu- 
lars, with a possible inclusion of blue elliptical systems. 



4 RESULTS 

The distribution of photometric redshifts in the cluster 
fields, after subtracting background galaxies, is shown in 
the left and middle panels in Fig. 2. The right-hand panel 
shows the redshift distribution of the background field. The 
location of the clusters is apparent for all clusters except 
C11600-I-41. We comment on this cluster in next subsection. 

As a measure of the richness of the clusters we give in 
Table 5 the number of cluster galaxies with Mg < — 17.7 
for the total areas (see Table 3), as well as inside a com- 
mon radius of 1 Mpc. We also list the density of background 
galaxies within the redshift range adopted for the differ- 
ent clusters. For the lower-redshift cluster, MS1008-1224, 
we also give numbers for Mg < — 16.2. 

In the left-hand panel of Fig. 3, we show the surface 
density of galaxies with Mg < — 17.7 for the four clus- 
ters as a function of radius. The horizontal solid line marks 
the surface density of the background field in the redshift 
range of C11600-I-41. Error bars and dashed lines represent 
Iff errors. The right-hand panel of Fig. 3 shows the projected 
fraction of early-type galaxies for the four clusters. The hor- 



izontal line marks the early-type fraction of the background 
field at « ~ 0.55, which is also similar to what we find at 
z ~ 0.31. All clusters have an early-type fraction clearly 
above the field value in the core region. In C10016+16 this 
fraction is ~85 per cent in the core. This decreases with ra- 
dius, but stays above 50 per cent for the whole area covered. 
The early-type fractions in C11601-h42 and MS1008-1224 is 
~60 per cent in the core, and show a decrease which ap- 
proaches the field value at large radius, even though the 
deviation from the field is still significant at the outermost 
points. For C11600-I-41, only the innermost point has a late 
type fraction deviating from the field. From both figures it 
is clear that we do not reach the field in any of the clusters, 
except C11600+41. 

In Fig. 4, we show the observed colour-magnitude (CM) 
diagram for all galaxies within the cluster redshift ranges for 
the three rich clusters, with different symbols representing 
early-type and late-type galaxies. At the redshifts of the 
clusters, the plotted colours approximately represent rest- 
frame B — V. From the figure, it is evident that the early- 
type galaxies populate a fairly narrow sequence in the CM 
diagram, which is clearly redder than the late-type galaxies. 
Also, the relative fraction of early-type galaxies decreases at 
fainter magnitudes. Note, however, that the spectral type is 
not independent of e.g., the rest-frame B — V colour, since 
colours are used when determining the photometric redshift, 
as well as spectral type. 



4.1 The poor cluster CI1600-|-41 

From the results above it is clear that C11600+41 is very 
poor, with only 36 ±24 cluster galaxies with m_R < 25. Fig. 
3 also indicates that if at all significant, the radial extent of 
the cluster is less than ~ 0.5 Mpc. At larger radii the surface 
density is marginally below the field value, while the late 
type fraction is consistent with the late type fraction of the 
field. 

The most significant indication of a cluster comes from 
the colour distribution, where the innermost point shows a 
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Figure 2. The left and the middle panels show the distribution of photometric redshifts in the cluster images after subtracting background 
galaxies. The right-hand panel shows the redshift distribution of the background field. Note the different scales on the y-axis. The peak 
in the background distribution at z ^ 1.4 is most likely caused by the lack of infrared photometry. All distributions above z ~ 1 
should therefore be viewed by caution. 




Radius (Mpc) 



Radius (Mpc) 



Figure 4. Left-hand panel: Projected surface density as function of radius for C10016+16 (filled circles), C11600+41 (squares), C11601-I-42 
(open circles) and MS1008— 1224 (triangles). The horizontal line shows the surface density of the background field in the redshift range 
of C11600-I-41. Right-hand panel; Early- type fraction as function of radius for the four clusters. The horizontal line shows the early-type 
fraction of the field at z = 0.54. The limiting magnitude is Mb = — 17.7. Error bars and the dashed lines represent Icr errors. 



clear excess of red galaxies (Fig. 3, right-hand panel). Also, 
the magnitude and colour of the central galaxy is similar 
to the brightest cluster galaxy (BCG) in the other clusters. 
Furthermore, a visual inspection of a three-colour image of 
C11600-I-41 shows this bright red central galaxy to be sur- 
rounded by a number of red and blue galaxies, clearly indi- 
cating the presence of a cluster. 

If we calculate the number of cluster galaxies in 
C11600+41 by the standard method of subtracting galax- 
ies in a blank field from the cluster image in a single band, 
we get a total of 9 ± 127 cluster galaxies to niR < 25. The 
large error is due to field-to-field variations, which dominate 



in the subtraction method. It is obvious that the subtraction 
method would not result in any cluster detection when using 
our blank field as a reference for the background counts. 

With only 36 selected cluster members within the cho- 
sen redshift range, the background counts dominate over the 
cluster galaxies, i.e. only 22 per cent of the galaxies are ex- 
pected to be cluster galaxies, as compared to 78, 66 and 80 
per cent, for C10016-M6, C11601-F42 and MS1008-1224. It 
is clear that this cluster is close to the detection limit. The 
large contamination in C11600-I-41 makes it impossible to 
determine either the internal properties of the cluster with 
any significance, or the cluster LF. Therefore, we do not in- 
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Figure 3. Obsereved colour-magnitude diagrams for C10016+16 

(top), C11601+42 (middle) and MS1008-1224 (bottom). Galaxies 
are separated into early-types and late-types. At the redshifts of 
the clusters, the observed colours approximately match rest- frame 
B — V colour. 



elude C11600-I-41 when comparing the LF and the Butcher- 
Oemler (BO) effect between the clusters. 



4.2 The cluster luminosity functions 

When calculating the LF we divide the galaxies within the 
cluster redshift range into magnitude bins with Am — 0.5 
to Mb = —17.7 for the two clusters at 2; ~ 0.55, and to 
Mb = -16.2 for MS1008-1224 at 2; ~ 0.31. K-corrections 
are determined for each galaxy using the best-fitting tem- 
plate from the photometric redshift calculations. For each 
magnitude bin we subtract background galaxies within the 
redshift range of the different clusters. 

The left-hand panels in Fig. 5 show the resulting LFs 
for the three clusters. As a first step we fit the total LFs to 
the usual Schechter function (Schechter 1976), 



$(M) 



10 



-0.4(a+l)M 



(2) 



where M* is the characteristic magnitude, representing the 
turnolf at the bright end of the LF profile, and a is the slope 
at the faint end of the LF. As an alternative to a, we also 



calculate the slope of a straight line fit to the five faintest 
bins in each LF, according to 

-1)M 



$/(M) oc 10-°-*("/^ 



(3) 



(Trentham 1998a; DFN02). The parameter a/ gives a bet- 
ter representation of the faint-end of the LF, since it is not 
affected by the coupling between M* and a in the Schechter 
function. In Table 6 we list M* , a and q/ for the clusters. 
Parameters are given both for the total cluster populations, 
and for the late-type population in each cluster. 

Fig. 5 shows that the shapes of the LFs of the three clus- 
ters differ substantially. To further understand these differ- 
ences, we plot the LFs divided into early-type and late-type 
galaxies separately in the right-hand panels of Fig. 5. For all 
three clusters, the early-type galaxies with Mb ^ — 18 have 
a Gaussian LF, peaking at Mg ~ — 20 (Table 6). We note 
that the faint early-type population in MS1008-1224 shows 
an increase below Mb = — 18. We return to this point in 
Section 5. 

The late-type galaxies are better fitted by Schechter 
functions. However, the relative normalisations of the early- 
type and late-type LFs vary between the different clusters. 
In Section 5 we discuss these differences further. 



4.3 The Butcher-Oemler effect 

We calculate the blue fraction, /b, using the definition in 
B084. For each cluster we select galaxies with My < — 20 
within a circular radius, -R30, containing 30 per cent of the 
cluster galaxies. Now, Jb is the fraction of these galEixies 
that are at least 0.2 mag bluer than the colour-magnitude 
relation for the cluster. To account for cluster galaxies out- 
side our field-of-view, we fit the cluster radial profile with a 
singular isothermal sphere when calculating R30. 

In order to compare results, we adopt in this section 
only the cosmology used by B084 (go = 0.1). The difference 
from B084 is that we here select cluster members with the 
photometric redshift method. This increases the accuracy 
since background subtraction is reduced. The blue fractions 
derived axe given in Table 7. 

For MS1008-1224, we find a value of /b that matches 
the straight-line fit of the blue fraction as a function of red- 
shift presented in B084. For C11601+42 and C10016-H16 the 
values of /b are located below this fit. The scatter in the 
fit is, however, large, and at least for C11601-f42 this devia- 
tion is hardly significant. C10016-I-16 is included in the B084 
sample, and they also note that this cluster is exceptional in 
the sense that it has a fraction of blue galax;ies that is more 
similar to local clusters than to clusters at « > 0.3. The 
blue fraction found by B084 is /b = 0.02 ± 0.07, which is 
consistent with our estimate. 



4.4 Brightest cluster galaxies 

In Table 7 we list the absolute B magnitude and rest-frame 
B — V colour for the brightest cluster galaxy (BCG) in 
the observed clusters, including the poor cluster C11600-I-41. 
Even though our sample is small, the results are consistent 
with the magnitude of the BCG being independent of the 
richness of the cluster, as shown by e.g. Sandage (1976) 
and Postman & Lauer (1995). In particular, we note that 
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Figure 5. Luminosity functions in the rcst-framc B-band for C10016+16, C11601+42 and MS1008-122'1. The left-hand panels show the 
total LFs, while in the right-hand panels the LPs are divided into early-type (open circles) and late-type (filled circles) galaxies. Error- 
bars include Poissonian uncertainties and field-to-field variance. For MS1008-1224 we show separate Schechter fits for the early-type and 
late-type populations. 



Table 6. Characteristic magnitude, Mg and faint-end slope, a derived from the Schechter function fit to the rest-frame B— band LPs 
for the total population and for late- type galaxies only, a/ is the slope of a straight-line fit to the five faintest-magnitude bins in each 
cluster. Limiting magnitude is = — 17.7. For MS1008-1224, we also give results to Mg = — 16.2. For the early-type population, 
is the peak magnitude and a'^ is the width of a Gaussian fit to Mb < — 18.0. 

Cluster All galaxies Late-types Early-types 

Mg a Of Mg a af a'^ 

C10016+16 -21.07 ±0.24 -0.55 ±0.10 -0.65 ±0.10 -21.75 ±0.50 -1.12 ±0.17 -1.17 ±0.18 -20.27 ±0.08 1.01 ±0.05 

C11601+42 -21.87 ±0.36 -1.25 ±0.10 -1.33 ±0.14 -21.80 ±0.53 -1.53 ±0.27 -1.65 ±0.19 -20.05 ±0.17 1.30 ±0.10 

MS1008-1224 -21.06 ±0.17 -0.72 ±0.26 -0.99 ±0.15 -21.05 ±0.58 -1.11 ±0.33 -1.61 ±0.24 -20.46 ±0.10 0.73 ±0.06 

MS1008-1224" -21.48 ±0.17 -1.46 ±0.09 -1.92 ±0.13 -22.00 ±0.37 -1.50 ±0.17 -1.89 ±0.16 

Note: 

a) Limiting magnitude Mg = — 16.2. 



the BCG of the extremely poor cluster C11600-h41 is only 5 DISCUSSION 
marginally fainter than the rest, demonstrating that even in 

this environment a luminous elliptical can be formed. g ^ Variations in the cluster luminosity function 

The different shapes of the LFs in our sample suggest that 
there is no universal form of the total cluster LF at « ^ 0.3 
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Table 7. The blue fraction, fg, as defined by B084, the absolute B magnitudes and rest- frame B — V colour 
for the brightest cluster galaxy in our cluster sample and the early-type fraction at two limiting magnitudes. 

Cluster 2 fs Mb B - V 

C10016+16 0.546 0.04 ± 0.02 -23.52 ±0.02 1.00 ±0.03 

C11600+41 0.540 - -23.36 ±0.03 0.92 ±0.05 

C11601+42 0.539 0.15 ±0.06 -23.50 ±0.04 0.93 ±0.07 

MS1008-1224 0.306 0.13 ±0.06 -23.65 ±0.05 1.15 ±0.07 



(Fig. 5). To quantify this we construct an average LF from 
the three individual best-fitting LFs. We then calculate the 

reduced chi-squarc when fitting each LF to the average LF. 
We also calculate the reduced chi-square when fitting the in- 
dividual LFs to a Gaussian+Schechter function. Results are 
shown in Table 8. The chi-squares are significantly higher 
for the fits to the average LF, supporting our claim for a 
non-universal LF. We also find that a single Schechter func- 
tion gives a poor representation of the cluster LF, which 
agrees with e.g. Driver et al. (1994) and Wilson et al. (1997). 
This is especially evident for MS1008-1224, where the best 
Schechter fit to Mb < — 17.7 yields a faint-end slope that 
is decreasing, while the data-points clearly show an increase 
in the number of faint galaxies with Mb > — 19. As in the 
case for Virgo, a better representation is given by the sum of 
a Gaussian and Schechter function. The reduced chi-squarc 
for for a Gaussian+Schechter function is x^/'^ = , 
while a single Schechter function yields x^/^* = 5.8. 

From the separate late-type and early-type LFs (right- 
hand panels of Fig. 5), it is clear that the relative abundance 
of the early-type and the late-type populations determines 
the over-all shape of the LF. In C11601+42 late-type galax- 
ies dominate at Mb ^ — 20. This results in a total LF 
that increases over the whole magnitude range, and has a 
steep faint-end slope. MS1008-1224 is similar to C11601+42 
in that it contains a numerous population of faint late-type 
galaxies, which dominates the total LF at faint magnitudes. 
The magnitude where this cross-over takes place is, however, 
slightly fainter in MS1008-1224 [Mb ~ - 19.5) compared 
to C11601+42 {Mb ~ - 20), which leads to a total LF 
where both an intermediate-magnitude Gaussian part, as 
well as a steep faint-end slope of late-type galaxies are dis- 
tinguishable. Note also that MS1008-1224 has a population 
of faint early-type galaxies, which adds to this. 

In C10016+16 there axe relatively few faint late-type 
galaxies, compared to the number of early-type galaxies at 
-21 ^ Mb ^ - 19. This leads to a total LF that has a 
Gaussian shape at intermediate magnitudes, and only at the 
faintest bin is there an indication of a rise. 

As discussed in DFN02, the differences between the 
shapes of the cluster LFs can be explained by the fact 
that clusters in the hierarchical clustering scenario at high 
redshift have a larger fraction of newly accreted galaxies 
with ongoing star formation. A general increase in the star- 
formation rate in field galaxies with redshift (Diaferio et al. 
2001), also contributes to an increase of star- forming galax- 
ies in clusters at high z. 

After accretion, and a possible period of enhanced star 
formation, the galaxies could have their gaseous envelopes 
removed by tides or ram pressure stripping, leading to a fad- 
ing over time-scales of Gyrs, as the remaining gas reservoir is 
exhausted ('strangulation') (Balogh et al. 2000; Diaferio et 
al. 2001). Such fading of the faint blue galaxies is discussed 



Table 8. Reduced chi-squares when comparing the LFs of 
C10016+16, C11601+42 and MS1008-1224 to the best-fitting 
Gaussian+Schechter function, as well as an LF which is an aver- 
age of the three individual LFs. 



Cluster 

C10016+16 
C11601+42 
MS1008-1224 



(individual fit) 
0.45 
0.46 
0.98 



(average fit) 
1.25 
1.49 
1.46 



by Wilson et al. (1997), who find that a fading of dwarf ir- 
regulars by ~ 3 mag can explain the difference between the 
LF in clusters z ~ 0.2 and the local Virgo cluster. Further, 
Conselice et al. (2001) show that spirals accreted at high 
z can be transformed into the faint dE population seen in 
local clusters by "galaxy harassment" (Moore et al. 1998). 

It is of obvious interest to compare the LFs of the differ- 
ent populations of these clusters with those of nearby clus- 
ters like Coma and Virgo. We then have the paradoxical sit- 
uation that studies of nearby cluster are limited to one, two 
or at most three colours, or alternatively to a morphological 
separation of different classes. A direct comparison with our 
spectral classification is therefore difficult. As a first step we 
therefore compare the total LF of our clusters with that of 
Coma at z = 0.02 (Trentham 1998a) and Virgo z = 0.003 
(Trentham & Hodgkin 2002). 

In Fig. 6 we plot the rest-frame S-band LF for 
C11601+42 (open circles), C10016+16 (filled circles), and 
MS1008-1224 (triangles) together with the LF for the rich 
Coma cluster (squares) and the poor Virgo cluster (dia- 
monds). For clarity, the LFs have arbitrary off-sets in the 
y-direction. As suggested by our discussion of the individual 
LFs, we represent the total LF as a sum of a Gaussian and a 
Schechter function. Because the bright and faint populations 
have distinctly different properties, we discuss them below 
separately. 



5.1.1 The bright population 

For Coma the Gaussian peaks at Mb = —19.3 ± 0.1, indi- 
cating a fading of the bright population by AM ~ 1 mag, 
compared to the high-redshift clusters (Table 6). The total 
LF of the Virgo cluster only has a marginal Gaussian peak. 
However, the population of morphologically classified ellip- 
tical galaxies, as well as giant spirals, has a clear Gaussian 
shape peaking at Mb ~ — 19 (Binggeli et al. 1988; Fergu- 
son & Sandage 1991), supporting a fading of this population 
similar to that suggested for Coma. 

These results are consistent with Small et al. (1997), 
who find that the characteristic magnitude My for early- 
type galaxies fade by AM ~ 0.7 mag between z = 0.54 
and z = 0. (We have here converted the results in Small 
et al. to our adopted cosmology.) Small et al. find a weaker 
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trend for bright late-type galaxies, which fade by AM<^ 0.4 
mag. An evolution is also found by Kodama & Bower (2001), 
who estimate that the bright blue galaxy population fades 
by AM ~ 1 mag between z ~ 0.4 and z = 0. De Propris 
et al. (1999) investigate the K band luminosity evolution of 
the bright galaxy population in clusters at 0.1 < 2; < 1, 
and find that K''{z) is consistent with passive luminosity 
evolution. We have compared our luminosity evolution with 
stellar synthesis models from Bruzual & Chariot (2003). We 
model the Gaussian population with an carly-typc popu- 
lation characterized by a single burst stellar population of 
age 8.0 Gyr and a Salpeter IMF. This combination results 
in colours that match the carly-typc spectral energy distri- 
bution in Coleman et al. (1980). For passive evolution, we 
find that this population should fade by ~ 0.7 mag since 
z ~ 0.55. This is slightly less than the ~ 1 mag evolution 
that our results suggests. This difi^erence will diminish if we 
either assume a Scalo IMF. or that the Gaussian population 
includes a fraction of galaxies with residual star formation. 



5.1.2 The faint population 

In DFN02 we argued that the different shapes of the LF 
at Mb ^ — 20 between Coma and C11601+42 can be un- 
derstood as a result of the dynamically younger age of the 
latter cluster. The steep blue end of the LF in this cluster 
should then consist of recently accreted field galaxies, which 
is supported by the fact that this part of the LF is almost 
exclusively made up of late- type galaxies. The rapidly in- 
creasing fraction of late-type galaxies in the outer parts of 
the cluster is consistent with this accretion scenario, and is 
also refiected in the rising blue fraction at large radii in this 
cluster, as shown in Section 5.4. A fading of the blue star 
forming galaxies in C1160H-42 by ~ 2 mag, would trans- 
form the LF of C11601-I-42 into a LF similar to Coma. 

The steepening of the faint-end in MS1008-1224 is 
clearly shown in Fig. 6. The slope is somewhat steeper than 
for Coma, and is shifted towards brighter magnitudes. From 
Fig. 5 wc see that the faiut-cnd of MS1008 1224 is made up 
of both blue and red gala^xies. The slope of the faint blue 
population is similar to the blue population in C11601-I-42, 
and a fading by ~ 1 mag of this population would result in a 
LF similar to the one in MS1008-1224. A subsequent fading 
by an additional magnitude will make both C11601-I-42 and 
MS1008-1224 similar to Coma. The fading of the late-type 
population between z = 0.5 and z = for these two clus- 
ters could therefore be described by a relation Am ~ — 4« 
mag. 

5.1.3 The red cluster C10016+16 

The LF of C10016-hl6 at 2 = 0.55 differs clearly from 
CI 160 1+42 at the same rcdshift. The LF of the former is 
similar to Coma in the plateau region, having a Gaussian 
shape, but shifted to brighter magnitudes by ~ 1 mag. From 
Fig. 5 we see that the relatively low abundance of late-type 
galaxies in C10016-I-16, as compared to C11601-I-42, is re- 
sponsible for the different shapes of the LF at intermediate 
magnitudes. 

Previous investigations of C10016-I-16 show that the 
bright galaxy population is dominated by red early-type 



galaxies. Koo (1981) found that this cluster has an unusu- 
ally high fraction of red galaxies compared to other high- 
redshift clusters, and suggested that the star formation must 
have ended a few Gyrs before the observed epoch. From 
the absence of blue galaxies and the small scatter in the 
colours of the red population. Small et al. (1995) conclude 
that the cluster is old, despite its high redshift. Of the ten 
clusters in the MORPHS sample (Small et al. 1997) with 
morphological classification from HST in the redshift range 
0.37 < z < 0.56, C10016+16 has the lowest fraction of 
spiral galEixies, fsp = 21 per cent, compared to a mean of 
fsp = 44±7 per cent for the remaining nine clusters. 

The low spiral fraction is, not unexpectedly, related to 
the low blue fraction found in C10016+16, which indicates 
that the cluster does not follow the general blueing described 
by the BO-effect. This may be explained by a higher dynam- 
ical age for C10016-I-16, which is supported by the high X-ray 
luminosity and velocity dispersion (Table 3), indicating that 
C10016-I-16 is more relaxed than C11601+42. 

C10016-I-16 has a very high fraction of post-starburst 
galaxies. These galaxies have no current star formation, but 
were forming stars ~l-2 Gyr before we observe them. Pog- 
gianti et al. (1999) found that 32±9 per cent of the galaxies 
brighter than My ~ —20.5 in C10016-f 16 are spectroscop- 
ically consistent with this type of galaxies. This suggests 
that a large fraction of the galaxies in C10016-I-16 has expe- 
rienced star formation since z ^ 0.8, but have subsequently 
faded to become red aX z = 0.54. Post-starburst galaxies 
are less frequent in local clusters (Poggianti et al. 1999), and 
although not actively forming stars, C10016-I-16 is therefore 
not similar to a local cluster placed at high redshift. This is 
also supported by the brighter early-type population in this 
cluster compared to local clusters. 

The difference between rich clusters at high and low 
redshift can qualitatively be explained within the hierarchi- 
cal clustering scenario. Kauffmann (1995) shows that high-z 
rich clusters assemble over a shorter time interval than \aw-z 
clusters of similar richness. In general, this naturally leads 
to a higher fraction of blue galaxies in high-z clusters. The 
properties of C10016-fl6 can then be understood if it formed 
during a short time interval at z 0.8, leaving a high frac- 
tion of post-starburst galaxies. This suggests that the cluster 
was blue at 2 ~ 0.8, while the red colours at 2 = 0.55 
indicates a very low infall of field galaxies during the last 
~ 1 Gyr before it is observed. 

It is possible that the special environment of C10016+16 
could have an influence on the blue fraction and the clus- 
ter LF. Connolly et al. (1996) and Hughes & Birkinshaw 
(1998) provide strong evidence that C10016-I-16 is part of a 
supercluster structure, with two associated clusters at simi- 
lar redshift with projected distances of 5 and 13 Mpc. The 
formation of this giant structure could have depleted the 
number of field galaxies surrounding C10016-I-16, and there- 
fore the accretion rate of late type galaxies. 

5.1.4 A universal cluster LF? 

To summarise our discussion, we find no evidence for a 
universal shape of the total cluster LF, in agreement with 
e.g., Binggeli et al. (1988) and Driver, Couch & Phillipps 
(1998). The LFs for the spectroscopic early- type and late- 
type galaxies separately have, however, similar shape in the 
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Figure 6. Luminosity function in the rest-frame B— band for 

C11601+42 at z = 0.54 (open circles), C10016+16 at z = 0.55 
(filled circles), MS1008-1224 at z = 0.31 (triangles). Coma at 
2 = 0.02 (squares) and Virgo at 2 = 0.003 (diamonds). Data for 
Coma and Virgo are taken from Trentham (1998a) and Trentham 
& Hodgkin (2002), respectively. The LFs are arbitrary off-set in 
the jy-direction. 



different clusters, but with varying relative strengths. This 
suggests an universality for type-specific LFs. This is analo- 
gous to the claim of type-specific LFs for different morpho- 
logical types by Biuggeli et al. (1988) and Andreon (1998). 
A comparison with Binggcli et al. (1988) is especially inter- 
esting, since they find that the bright part of the LF mainly 
consists of elliptical galaxies with a Gaussian LF, and that 
the steep faint-end consists of irregulars and dEs. This sug- 
gests, as expected, that the early-type population mainly 
consists of ellipticals, while many faint late-type galaxies 
are irregulars. To explain the different shapes of the high-« 
cluster LFs, we argue that these reflect different dynamical 
states of the clusters. As the clusters get dynamically older, 
we expect the late-type population to fade relative to the 
early-type population, and LFs to become more similar. 



5.2 The dwarf population in MS1008-1224 

The faint limit in absolute magnitude we reach for MS1008- 
1224 at 2 = 0.31 allows us to study the dwarf popula- 
tion ^ 1.5 mag deeper in this cluster compared to the 
clusters at z ~ 0.55. Studies by Trentham (1998c) and 
Boyce et al. (2001) have shown that dwarf irregular (dirr) 
galaxies have B — R ~ 0.9, while dEs have colours in a 
broader range l.Z <^ B — R <^ 2.0. To compare the colours 
of the dwarf population in MS1008-1224 with other nearby 
clusters, with more limited colour information, we therefore 



plot in Fig. 7 the rest-frame B — R colours of galaxies with 
-17.7 < Mb < - 16.2. 

In t\ieB-R histogram of MS1008-1224 in Fig. 7 there 
are two peaks, suggesting that there indeed are two distinct 
populations representing dirrs and dEs. In the figure we also 
show the distribution ol B — R colours for the faint galaxies, 
which on the basis of the full BVRI photometry we have 
classified as early-type and late-type galaxies, respectively. 
The fact that this classification closely follows the two-colour 
division into the two peaks shows that wo reliably can use 
the B — R colour to broadly distinguish between early-type 
and late-type galaxies. 

From the B — R colours Trentham (1998c) finds that 
all dwarf galaxies with -18.9 < Mr < - 16.9 in Abell 
665 ai z = 0.18 arc consistent with being dEs (called dSphs 
by Trentham). In Abell 963 a,t z = 0.21 the majority of 
the dwarf galaxies with —19.2 < Mr < — 18.2 also 
have colours consistent with dEs. There is, however, in this 
cluster also a population with colours intermediate between 
dEs and dlrrs, which Trentham proposes may represent a 
transitional stage between these types. 

In a recent study Boyce et al. (2001) investigate the 
_B — _R colour of galaxies with —19 < Mr < —16.5 in Abell 
868 aX z = 0.154. They find a distribution in colour with 
two peaks, showing the presence of a dominant population of 
dEs, but also a population of dlrrs. For the faintest galaxies 
(—17.5 < Mr < — 16.5) there are similar numbers of dEs 
and dlrrs. 

The main difference between MS1008-1224 and lower- 
redshift clusters is the relative fraction of dlrrs to dEs. To 
make a quantitative comparison between the faint popula- 
tions, we use the B — R index to calculate the fraction of 
dlrrs for MS 1008-1224, AbeU 963. Aboil 665 and Abell 868. 
We assume that galaxies with rest-frame B — R < 1.2 are 
dlrrs, which is the criterion used by Boyce et al. (2001). 

A complication is that MS1008-1224 is not observed 
in standard Johnson-Cousins filters, but in Bessel filters. 
To correct for this we use galaxy templates from Coleman 
et al. (1980) and filter transmission curves to calculate the 
transformation between the standard B — R colours and the 
B — R colours for the Bcssel filters. We then find that the 
division between dlrrs and dEs B — R = 1.2 in standard 
Johnson-Cousins filters corresponds to B — R = 1.34 in the 
Bessel filters. 

Using this criterion, we find that the fraction of faint 
galaxies that belong to the dIrr population is 0.72 ± 0.03 
for MS1008-1224, to be compared with 0.07 ±0.03 for Abell 
963, 0.01 ±0.01 for Abell 665 and 0.22 ±0.02 for Abell 868. 
A further problem when comparing these numbers is that 
the different investigations use different magnitude intervals 
for defining the dwarf populations. Irrespective of this, how- 
ever, we reach the important conclusion that MS1008-1224 
is dominated by dlrrs, while the lower redshift clusters are 
dominated by dEs. 

The fraction of dlrrs in Coma, classified by the 
same colour criterion as above, for galaxies with 
-15.3 < Mr < -13.3 is 0.11 ±0.01 (calculated from fig. 
7 in Trentham 1998a). The fraction of morphological dlrrs 
in Virgo ( — 17.8 < Mr < — 15.8), calculated from fig. 
6 in Trentham & Hodgkin (2002), is 0.25±0.07 The results 
emphasize the conclusion that dEs dominate the faint pop- 
ulation in nearby clusters, while our study shows that the 
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Figure 7. The rest- frame B — R colour distribution for galaxies in 
MS1008-1224 with -17.7 < Mb < - 16.2. The solid line shows 
the total number of galaxies, while the dotted and dashed lines 

show the distribution of carly-typc and latc-type populations, as 
classified by multi-band photometry, separately. 



opposite is true for MS1008-1224 at z = 0.31. It would ob- 
viously be of great interest to study more clusters to these 
limits to see if this is a general property. 

Galaxy harassment could here bo a mechanism that 
transforms the faint late-type population, dominating 
medium redshift clusters such as MS1008-1224, into dEs, as 
shown by Conselice et al. (2001). The similarity between the 
faint-end slopes of the late-type population (a/ = — 1.84± 
0.12) and the carly-typc population (a/ — — 1.69 ± 0.14) 
at faint magnitudes {Mb > — 19) is consistent with this 
picture. Based on these observations alone, we cannot, how- 
ever, rule out that slower mechanisms such as 'strangulation' 
(Balogh et al. 2000), are at least partly responsible for the 
luminosity and colour evolution. 



5.3 Dwarf-to-giant ratio 

Previous studies have shown that the ratio of faint to bright 
galaxies increases at large radii, where the surface density 
decreases. Driver et al. (1998) find that five out of seven 
clusters at z ~ 0.15 show an increasing dwarf-to-giant 
ratio with radius. A K-hand study of AC 118 at 2 = 0.3 
by Andreon (2001) shows the same trend. In DFN02 we 
found that C11601-I-42 has a steeper faint-end slope of the 
LF in the outer part of the cluster, which is equivalent to 
an increasing dwarf-to-giant ratio. 

Here, we calculate the dwarf-to-giant ratio by defining 
galaxies with —19.5 < Mb < — 17.7 as dwarfs, and 
brighter galaxies as giants. The separation between dwarfs 
and giants is chosen to match the cross-over in the clus- 
ter LPs between early-type and late- type galaxies (Fig. 5). 
This division is somewhat brighter than used in other papers 
(e.g.. Driver et al. 1998), and is here adopted in order to get 
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Figure 8. Dwarf-to-giant ratio as a function of bright galaxy 
surface density for C10016-I-16 (circles), C11601-I-42 (squares) and 
MS1008— 1224 (triangles). Filled symbols represent the inner re- 
gion of the cluster (R < 0.5 Mpc), while open symbols represent 
the outer region. The horizontal lines represent the field value 
(solid line) and lu-errors (dotted lines). 



sufficient statistics in the faint bin. We further divide each 
cluster into a core region with projected radius < 0.5 Mpc, 
and an outer region between 0.5 and 1.1-1.5 Mpc, depending 

on cluster. 

In Fig. 8 we plot the resulting dwarf-to-giant ratios as 
function of surface density of the giant galaxies for these two 
regions in all clusters. C10016-fl6 is represented by circles, 
C11601+42 by squares and MS1008-1224 by triangles. Filled 
symbols show results for the core area, while open symbols 
represent results from the outer area. The horizontal lines 
represent the field value and the corresponding la-errors. 

Pig. 8 shows that C10016-M6 and C11601-F42 have a 
clear trend of an increasing dwarf-to-giant ratio when go- 
ing from the inner high-density region, to the outer lower- 
density region, which is consistent with the trend found by 
Phillipps et al. (1998) and with Andreon (2001). As an ex- 
planation of this these authors suggest that the high-density 
environment in cores of clusters is hostile to the dwarf galax- 
ies, while the outer, low-density regions do not affect this 
population, which therefore have a dwarf-to-giant ratio and 
faint-end slope similar to that of the field. 

For the low redshift cluster MS 1008-1224, there is no 
obvious trend. We note, however, that the error in the dwarf- 
to-giant ratio is substantial. 



5.4 The Butcher— Oemler effect as function of 
radius and limiting magnitude 

The BO effect is of major observational importance as a 
probe of the cluster galaxy population. Although much of 
the discussion below to a large extent is a consequence of the 
early-type and late-type LPs in Pig. 5 and their dependence 
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Figure 9. The blue fraction, fs, as a function of radius for 
C10016+16 (filled circles), C11601+42 (open circles) and MS1008- 
1224 (triangles). R30 is the fiducial value used by B084 for cal- 
culating fs ■ 
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Figure 10. The fraction of blue galaxies, /s, as a function of lim- 
iting magnitude for C10016+16 (filled circles), C11601-I-42 (open 
circles) and MS 1008-1224 (triangles). 



on radius (e.g.. Fig. 3), we will for this reason illustrate the 
consequences of this with two important implications for the 
BO effect. 

Already in B084, and more recently in DFN02, it was 
pointed out that the blue fraction Jb increases at larger 
radii. This can be understood as a result of the fact that 
carly-typcs are in general more centrally concentrated than 
late-types. In Fig. 9 we show the dependence of fs on radius, 
where the radius is given in fractions of i?3o (as defined in 
Section 4.3). Also in this respect 010016+16 deviates from 
the other clusters. Inside of 1.257i3o, /s decreases rapidly, 
while outside this radius /s stays almost constant, at a value 
well below that of the other clusters. 011601+42 on the other 
hand shows a steep radial dependence. The higher blue frac- 
tion in 011601+42 is consistent with the higher fraction of 
late type galaxies in this cluster (Fig. 3). MS 1008-1224 has 
a trend in between the two higher rcdsliift clusters. From 
a very low fs in the core region, there is a rapid increase 
to ~ 1.25 R30, followed by a slower increase in the outer 
regions. The behaviour of /s in these clusters reflects, as 
expected, the radial dependence of the early-type fraction 
shown in the right panel of Fig. 3. 

In DFN02 we showed that fs in 011601+42 depends 
strongly of the limiting magnitude. This dependence is a di- 
rect consequence of the different shapes of the early- type and 
late-type LFs in this cluster. At fainter magnitudes the to- 
tal LF becomes more dominated by late-type galaxies, which 
leads to a higher /u. In Fig. 10 we plot the dependence of 
/b on the limiting magnitude for the three clusters. As ex- 
pected, the increase of /s is much weaker in C10016+16 than 
in 011601+42, which is consistent with the LF for this clus- 
ter being dominated by early-type galaxies to fainter mag- 
nitudes (Fig. 5). MS1008-1224 again shows an intermediate 



behaviour. The strong increase in /s at faint magnitudes 
in MS 1008-1224 is caused by the dominating population of 
blue dwarf galaxies. 



6 SUMMARY 

This work represents one of the deepest studies of the LF for 

intermediate-redshift clusters, several magnitudes beyond 
L* . In addition, we have demonstrated the usefulness of pho- 
tometric redshifts for selecting clusters members, as well as 
for population studies of distant clusters of galaxies. As a 
result, while most previous determinations of distant clus- 
ter LFs have treated only the total LF, we have in this study 
been able to separate the early-type and late-type popula- 
tions. This has been possible, despite the moderate telescope 
size due exactly to the use of photometric redshifts. Com- 
pared to the usual subtraction method, we have been gaining 
a large factor in respect to the background contamination. 
Our main conclusions from this study are: 

• There is no universal shape of the total cluster LF at 
z ^ 0.3. 

• The early- type population has a Gaussian LF, while the 
late-type population is well fitted by a Schechter function. 

This suggests that the LFs for different spectral populations 
could be universal, while the total LF depends on the rela- 
tive abundance of these populations. 

• The evolution of the late-type galaxies is consistent with 
a fading by ~ 2 magnitudes between z ~ 0.55 and z = 0, 
while the early-type population fades by ~ 1 mag. This 
scenario suggests that the total LFs of the high- 2; clusters 
become more similar to local LFs as the clusters get dynam- 
ically older. 
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• The red cluster C10016+16 is an atypical liigli-2 cluster 
that resembles local rich clusters in many aspects, indicating 
an old dynamical age despite its redshift. It does, however, 
contain a large fraction of post-starburst galaxies, suggesting 
that star formation was more intense at 2 J; 0.8, and that 
the infall was very low during the last ~ 1 Gyr before the 
cluster is observed. 

• In MS1008-1224 at z = 0.31, we find that dirrs domi- 
nate over dEs, opposite to what is found in nearby clusters. 
If this is confirmed for more clusters at this redshift it im- 
plies a dramatic evolution of the dwarf population. 

• The relation between dwarf-to-giant ratio and surface 
density indicates that high-density regions arc hostile to 
dwarfs, consistent with a destruction or fading of this pop- 
ulation by galaxy harassment. 

• We find that the blue fraction, /b, as defined by B084 
varies with radius and limiting magnitude. This is a direct 
consequence of the radial gradient of the late-type galaxies 
and the relative normalisation of the late-type and early- 
type LFs. 

There arc several natural next steps to this study. The 
sample needs to be expanded, both in terms of redshift 
and cluster properties. As the example of CI 0016-1-16 has 
shown, even at the same redshift the the population may 
vary greatly, depending on the dynamical state, and possi- 
bly on the environment. A better understanding of this is 
of obvious importance for the understanding of the cluster- 
formation process, especially at redshifts higher than we 
have probed. 

Spectroscopic studies at redshifts similar to the clusters 
in this study may today marginally be carried out with 8-10 
m class telescopes. This would of course give more detailed 
information about the dynamical properties of the clusters, 
as well as a general check of the reliability of the photometric 
redshift method for clusters. At higher redshifts it is, how- 
ever, virtually impossible to probe the population beyond L, 
with spectroscopy, even with these telescopes. For this type 
of investigations one therefore has to rely on photometric 
redshifts, in this case extended to include the near-infrared 
bands. The study in this paper, as well as DFN02, hopefully 
represents some first steps in this direction. 

Also at low redshifts there axe several open questions, 
especially connected to the dwarf population. The dramatic 
difference in the dIrr-to-dSph ratio seen between MS1008- 
1224 and the local clusters needs to be confirmed with more 
medium-redshift clusters. 

Finally, our study has mainly given information about 
the broad band colours of the cluster members. A high- 
resolution study with HST of the morphological properties 
of these galaxies would be of prime interest, as was partially 
demonstrated for CI 160H-42 in DFN02, and in more de- 
tail for the bright population by e.g. the MORPHS sample 
(Small et al. 1997). We hope in the future to pursue some 
of these questions. 
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